Discrete breathers in a nonlinear electric line: 
Modeling, Computation and Experiment 
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We study experimentally and numerically the existence and stability properties of discrete 
breathers in a periodic nonlinear electric line. The electric line is composed of single cell nodes, 
containing a varactor diode and an inductor, coupled together in a periodic ring configuration 
through inductors and driven uniformly by a harmonic external voltage source. A simple model for 
each cell is proposed by using a nonlinear form for the varactor characteristics through the current 
and capacitance dependence on the voltage. For an electrical line composed of 32 elements, we find 
the regions, in driver voltage and frequency, where n-peaked breather solutions exist and character- 
ize their stability. The results are compared to experimental measurements with good quantitative 
agreement. We also examine the spontaneous formation of n-peaked breathers through modula- 
tional instability of the homogeneous steady state. The competition between different discrete 
breathers seeded by the modulational instability eventually leads to stationary n-peaked solutions 
whose precise locations is seen to sensitively depend on the initial conditions. 

PACS numbers: 05.45.Yv, 63.20.Pw, 63.20.Ry 
Keywords: Localized modes; Solitons; Discrete Breathers 



I. INTRODUCTION 

Nonlinear physics of discrete systems has witnessed 
enormous development in the past years. In particular, 
a great deal of interest has been paid to the existence 
and properties of intrinsic localized modes (ILMs), or 
discrete breathers, which result from the combination of 
nonlincarity and spatial discreteness. These spatially lo- 
calized states have been observed in a wide variety of 
different systems [l], . They were originally suggested 
as excitations of anharmonic nonlinear lattices [3|, but 
the rigorous proof of their persistence under general con- 
ditions 3 led to their investigation in a diverse host 
of applications. These include, among others, antifcr- 
romagnets p|, charge-transfer solids Q, photonic crys- 
tals 0], superconducting Josephson junctions Q, mi- 
cromechanical cantilever arrays Q , granular crystals (lOj 



and biopolymers llj. More recently, the direct ma- 



nipulation and control of such states has been enabled 
through suitable experimental techniques fl2j . 

Despite the tremendous strides made in this field, the 
relevant literature, nevertheless, often appears to be quite 
sharply divided between theory and experiment. Fre- 
quently, experimental studies do not capture the dynam- 
ics in enough detail to facilitate an exact comparison with 
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theoretical studies. At other times, the theoretical mod- 
els are not refined enough (or lack the inclusion of non- 
trivial experimental factors some of which may be diffi- 
cult to quantify precisely) to make quantitative contact 
with the experimental results. Our system — a macro- 
scopic electrical lattice in which solitons have a time- 
honored history [l3|, [lij — is, arguably, ideally suited for 
this kind of cross-comparison: the lattice dynamics can 
be measured fully in space and time, and the physical 
properties of individual unit cells of the lattice can be 
characterized in enough detail to allow for the construc- 
tion of effective models. 

In this paper, we present a detailed study of discrete 
breathers in an electric lattice in which ILMs have been 
experimentally observed (l5l4l7| . We propose a theoret- 
ical model which allows us to systematically study their 
existence, stability and properties, and to compare our 
numerical findings with experimental results. We demon- 
strate good agreement not only at a qualitative but also 
at a quantitative level between theory and experiment. 
The presentation of our results is organized as follows. In 
the next section we study the characteristics (intensity 
and capacitance curves versus voltage) of the varactor, 
the nonlinear circuit clement, in order to develop the rel- 
evant model for the electrical unit cell. The results for 
the single cell are validated through the comparison of 
its resonance curves for different driving strengths. We 
also derive the equations describing the entire electrical 
line. In Sec. IIIII we study the existence and stability 
properties of n-peaked breathers for n = 1,2, 3 in the 
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FIG. 1: Schematic circuit diagram of the electrical transmis- 
sion line. 



driving frequency and voltage parameter space. The nu- 
merical results are compared to the experimental data 
with good quantitative agreement. We also briefly study 
the spontaneous formation of n-peaked breathers from 
the modulational instability of the homogeneous steady 
state. We observe that the location and number of the 
final peaks depends sensitively on the initial conditions. 
Finally, in Sec. IIV1 we conclude our manuscript and offer 
some suggestions for possible avenues of further research. 




FIG. 2: Experimental data (*) and numerical approximation 
thereof (continuous line) corresponding to C(V) and Id(V) 
(linear and semi-log plots) for the nonlinear varactor. 



II. THEORETICAL SETUP 

Our system consists on an electric line as represented 
in Fig. [TJ This line can be considered as a set of single 
cells, each one composed of a varactor diode (NTE 618) 
and an inductor Li = 330 /iH, coupled through inductors 
Li = 680 fiH. Each unit cell or node is driven via a resis- 
tor, R = 10 kCl, by a sinusoidal voltage source V(t) with 
amplitude Vd and frequency /. In experiments a set of 32 
elements have been used, with a periodic ring structure 
(the last element is connected to the first one), and mea- 
surements of voltages V n have been recorded. Related to 
the voltage source, we have considered amplitudes from 
Vd = 1 V to Vd = 5 V and frequencies from / = 200 kHz 
to / = 600 kHz. 

In order to propose a set of equations to character- 
ize the electrical line, we have used circuit theory and 
Kirchhoff's rules; the main challenge has been to de- 
scribe appropriately each element. In general, resistance 
and inductors are inherently imperfect impedance com- 
ponents, i.e., they have series and parallel, reactive, ca- 
pacitive and resistive elements. Moreover, due to the 
commercial nature of the elements, manufactured compo- 
nents are subject to tolerance intervals, and the resultant 
small spatial inhomogeneity introduces some additional 
uncertainty. We have quantified the spatial inhomogene- 
ity by separately measuring all lattice components. The 
diode capacitance was found to vary by 0.3 percent (stan- 
dard deviation), whereas the inductors both exhibited a 
0.5 percent variation. Additional factors that may con- 
tribute slightly to inhomogeneities are wire inductances 
as well as load and contact resistances. The var actors 
(diodes) we use are typically intended for AM receiver 
electronics and tuning applications. As described later, 



we characterize this lattice element in more detail, since 
it is the source of the nonlinearity in the lattice. 

As a guiding principle, we are aiming to construct a 
model which is as simple as possible, with a limited num- 
ber of parameters whose values are experimentally sup- 
ported, but one which is still able to reproduce the main 
phenomenon, namely nonlinear localization and the for- 
mation of discrete breathers. With this balance in mind, 
we proceed as follows. 

In the range of frequencies studied it is a good approx- 
imation to describe the load resistance as a simple resis- 
tor, neglecting any capacitive or inductive contribution. 
Also, we have performed experimental measures of the 
varactor characteristics. This experimental data shows 
that it can be modeled as a nonlinear resistance in par- 
allel with a nonlinear capacitance, where the nonlinear 
current Id(V) is given by 

I D (V) = -I a exv(-pV), (!) 

where /3 = 38.8 V" 1 and I a = 1.25 x 10~ 14 A (we consider 
negative voltage when the varactor is in direct polariza- 
tion), and its capacitance as 

(C V + C W (V') + C(V) 2 if V < v c , 
C(V) = { (2) 
{ C e~ aV if V>V C , 

where V = (V - V c ), C = 788 pF, a = 0.456 V" 1 , 
C v = Co exp(— aV c ), C w = —aC v (the capacitance and 
its first derivative are continuous in V = V c ), C = 100 nF 
and V c = —0.28 V. In Fig. [2] we present the experimental 
data and their corresponding numerical approximations 
for Id(V) and C(V), where a good agreement between 
the two can be observed. 
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With respect to the inductors, in the range of frequen- 
cies considered, capacitive effects are negligible, but they 
possess a small dc ohmic resistance which is around 251. 
The inductors and the varactor are a source of damp- 
ing in the ac regime, and these contributions must be 
taken into account. However, we have no manufacturer 
data related to dissipation parameters, and it is difficult 
to measure them experimentally. In order to introduce 
these effects, we will model dissipation phcnomcnolog- 
ically by means of a global term given by a resistance 
Ri, which appears in each unit cell in parallel with L 2 ; 
to determine its value, we have studied experimentally a 
single element as shown in Fig. [3] In this way, we will 
consider the inductors themselves as ideal elements. 
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FIG. 3: Single cell element model. 



Using basic circuit theory, the single element is de- 
scribed by the equations: 



dv 1 cos(57t) 

dr c(v) |_ RCqluq 

— 7=r ( -5- + d ) v + (y-io) 

ujqCq \Ri R 



dy ( L 2 . 



(3) 



(4) 



where dimcnsionless variables have been used: r = uiot, 
Id = Id/(woCoV,i), v — V/Vd, the dimcnsionless volt- 
age at point A, c(v) — C(V)/Cq, = oj/uo, and 
ujq = 1 / \JLiCq\ y represents the normalized current 
through the inductors. 

We can generate theoretical nonlinear resonance curves 
and, comparing with experimental data, select the opti- 
mal dissipation parameter value Ri . Results are summa- 
rized in Table HI and the comparison between theoreti- 
cal and experimental data is shown in Fig. 2) Also, we 
consider a small frequency shift of 18 kHz in numerical 
simulations to match the resonance curves. This effect 
may originate from some small capacitive and/or induc- 
tive contributions that we have not previously taken into 
account. 

Thus, with a consistent set of parameter values, and 
using again elementary circuit theory, we describe the full 
electric line of N coupled single cell elements by the fol- 
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TABLE I: Values of the resistance Ri corresponding to differ- 
ent voltage amplitudes for the driving source Vd- 
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FIG. 4: (Color online). Nonlinear resonance curves corre- 
sponding to a single element. Dots (orange [grey]) correspond 
to experimental data while the continuous and dashed black 
lines correspond, respectively, to stable and unstable numer- 
ical solutions. From bottom to top: Vd = 1, 2, 3, 4 and 5 V. 



lowing system of coupled ordinary differential equations: 



dv n . cos(fir) 

c{v n )— = y n -iD{v n ) + -— 

dr RCqojq 



1 1 \ v. 



Ri R) wqCo 



dy n L 2 
dr Li 



(l'n+1 + V n -l - 2v n ) - V n , (5) 



where all magnitudes are in dimcnsionless units. Within 
this nonlinear dynamical lattice model, our waveforms of 
interest, namely the discrete breathers, are calculated as 
fixed points of the map 



(G) 



where T = 1// is the temporal period of the breather. In 
order to study the linear stability of discrete breathers, 
we introduce a small perturbation {^ n ,T] n ) to a given so- 
lution (v nQ ,y n0 ) of Eq. © according to v n = v n0 + {;„, 
Dn — ynO + Vn- Then, the equations satisfied to first order 
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by (£n,Vn) are: 



c{v n0 ) 



dr 



T] n - T(v n0 ;t)£, n 



dr) n 



L2 

dr L\ 
with r(u n o;t) being 

diD(v n o) 



(£n+i + £„-i - 2£ n ) - £ n (7) 
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To identify the orbital stability of the relevant solutions, 
a Floquet analysis can be performed. Then, the stability 
properties are given by the spectrum of the Floquet oper- 
ator M. (whose matrix representation is the monodromy) 
defined as: 
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The AN x AN monodromy eigenvalues A are called the 
Floquet multipliers. If the breather is stable, all the eigen- 
values lie inside the unit circle. 



III. NUMERICAL COMPUTATIONS & 
EXPERIMENTAL RESULTS 
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Using Eq. ([5]), we have generated n-peak ILMs and 
determined numerically their stability islands in (V^, /) 
parameter space. We have found that there exist over- 
lapping regions where two, or several, of these n-peak 
configurations exist and are stable. Therefore, the long 
term dynamics in these regions is chiefly dependent on 
initial conditions. However, determining precisely the 
basins of attraction of each configuration is not possible 
because of the high dimensionality of the problem. 

Figure [5] shows the existence and stability diagrams 
for 1-, 2- and 3-peak breathers. The values of Ri are ob- 
tained through cubic interpolation from Table HI Notice 
that the lattice allows different configurations of solutions 
of n-peak breathers for n > 1 corresponding to the peaks 
centered at different locations, but determining precisely 
the diagrams corresponding to all possible configurations 
is not possible because of the size of the lattice. There- 
fore, we have studied only breathers with peaks as far 
apart as possible. 

Remarkably, breathers are generally robust for Vd < 
3 V. Above that critical value, the considered solutions 



FIG. 5: (Color online). Existence islands for n-peak 
breathers. The white line denotes the threshold where the 
homogeneous state becomes unstable (above the line). Black 
areas correspond to stable solutions and orange (grey) areas to 
unstable solutions, (a) one-peak breather (the homogeneous 
state threshold is located at the top border of the areas), (b) 
a family of two-peak breather and (c) a family of three-peak 
breathers. 



may become unstable in some "instability windows" (see 
orange [grey] regions). Those instabilities, which are of 
exponential kind, typically lead an onsitc breather pro- 
file to deform into a stable inter-site breather waveform. 
We have analyzed in more detail the dependence of those 
instabilities for I-peak breathers. The study for higher 
peaked structure is cumbersome due to the increasing 
number of inter-site structures that might arise. Fig- 
ure [S] shows an example of 1-peaked on-site and inter-site 
breathers together with their Floquet spectrum. 

A detailed analysis shows that for 3 < Vd < 5 V, inter- 
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FIG. 6: (Color online). Breather profiles (left column) and 
corresponding Floquet multiplier spectra (right column) for 
t = 0, V d = 4 V and / = 243.6 kHz. The top (bottom) 
panels correspond to the unstable (stable) one-site (inter-site) 
breather. 



site breathers are always stable, whereas on-site ones may 
be unstable (i.e., there is no stability exchange as it oc- 
curs in Klein-Gordon lattices). However, for < 2 V, 
inter-site breathers are always unstable whereas on-site 
breathers are stable. Finally, for 2 V< Vd < 3 V, inter- 
site breathers experience instability windows whereas 
their on-site counterparts are stable. Figures [THE illus- 
trate a typical set of relevant results for the cases of 
V d = 1.5 V, V d = 2.5 V and V d = 4 V which summa- 
rize the different possible regimes as Vd is varied. The 
unity crossings indicated by the red line (wherever rele- 
vant) in the right panel of each of the figures mark the 
stability changes of the pertinent solutions. 

We have also performed an analysis of the stability 
of 1-peak breathers for larger lattices (N = 101) and 
observe that the existence and stability range for on-site 
and inter-site breathers are not significantly altered. 

To corroborate the numerical picture, we have also 
performed an analogous experimental analysis. To deter- 
mine experimentally the stability islands, we have chosen 
a particular voltage Vd and, starting at low frequencies 
(homogeneous state), and increasing / adiabatically, we 
can reach different regions. Going up and down adiabat- 
ically we are able to examine regions of coexistence be- 
tween different n-peak breathers. In particular, regions 
corresponding to 1-pcak, 2-peak and 3-peak breathers 
have been detected experimentally. The 3-peak region is 
bounded from above by the 4- and 5-peak areas, which 
wc did not track. We have determined the boundaries by 
doing up-sweeps and down-sweeps, with significant hys- 
teresis phenomena. 

It should be mentioned that for the interior of the 1- 
peak region, we demonstrated experimentally that the 
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FIG. 7: (Color online). Dependence of the Floquet multi- 
pliers on the frequency for N — 32 cells and Vd = 1.5 V. 
The top (Bottom) row of panels corresponds to the on-site 
(inter-site) breather. The left and right column of panels de- 
pict, respectively, the argument and magnitude of the Floquet 
multipliers. The horizontal (red) line in the spectra depicts 
the stability threshold. 
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FIG. 8: (Color online). Same as Fig. but for V d = 2.5 V. 



discrete breather /ILM can be centered at any node in the 
lattice and survive there. This verification implies that 
the experimental lattice — despite its inherent compo- 
nent variability — does display a sufficient degree of spa- 
tial homogeneity for the basic localization phenomenon 
to be considered (discrete) translationally invariant. In 



6 



on-site breather 




f (kHz) f (kHz) 

FIG. 9: (Color online). Same as Fig. but for V d = 4 V. 
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FIG. 10: The ILM can exist at any site of the lattice; here f 
= 253 kHz and V d = 4V. We demonstrate that the ILM can 
be made to jump to the neighboring site upon the temporary 
creation of an impurity there. The left-most profile (centered 
at site 16) turns into the dotted (red) trace when an external 
inductor is placed in the direct vicinity of the L2 inductor 
at site 17. Upon removal of this external inductor, the ILM 
persists at site 17 (solid black trace). This process can be 
repeated to move the ILM further to the right, as shown. 



practice, we employed an impurity in the form of an ex- 
ternal inductor physically touching a L2 lattice inductor 
to make the ILM hop to that impurity site; see Fig. [TO] 
Upon removing the impurity, the ILM would then persist 
at that site. We believe that this technique could prove 
extremely valuable towards the guidance and manipula- 
tion (essentially, at will) of the ILMs in this system. 

A comparison between theoretical and experimental 
data of 1-peak and 3-peak breathers is shown in Fig. [TT] 
The precise width of any stability region is hard to match, 
because there exist both regions of different peak num- 
bers and regions of different families with the same num- 
ber of peaks overlap, and thus their competition prevents 
an absolutely definitive picture. Which one "wins" out 




FIG. 11: (Color online). Existence regions of 1-peaked (top 
panel) and a family of 3-peaked (bottom panel) breathers ob- 
tained numerically (black areas correspond to stable solutions 
and orange (grey) one to unstable solutions) as compared to 
experimental data identifying the range of observations of the 
corresponding type of states (circles). The theoretical data 
are displaced by a +7 kHz frequency offset. 



appears to depend sensitively on small lattice impurities 
present in our (commercial) experimental elements. For 
instance, in the experiment it looks as if the 2-peak region 
is squeezed in favor of the 1-peak region. This might be 
the reason why the experimental 1-pcak region is slightly 
wider at higher driver amplitudes than it appears in the 
corresponding theoretical predictions. Nonetheless, the 
comparison between experimental and theoretical exis- 
tence regions depicted in Fig. [TT] shows generally good 
qualitative (and even quantitative) agreement in the con- 
text of the proposed model. 

More detailed experimental results are shown in Fig.[T2l 
where we depict peak profiles at two different driver volt- 
ages. The profiles were taken at the times of largest peak 
voltage amplitude and lowest peak voltage amplitude. 
For Vd = 2V (left column of panels) we see that, as the 
frequency is raised from below, we cross from the 1-peak 
region through the 2-pcak region and into the 3-pcak 
region. For the Vd = 4 V case (right column of pan- 
els) the same sequence can be observed when scanning 
in one frequency direction. In order to illustrate both 
the hysteresis and the overlap between n-peak regions, 
we depict in the figure (panels (e) and (f)) a situation 
where the 2-peak solution occurs at a higher frequency 
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FIG. 12: Experimental (•) and numerical (o) 1-, 2- and 3-peak 
breather profiles for different driver voltage Vd and frequen- 
cies values. The left (right) column of panels corresponds 
to a driver strength Vd = 2 V (Vd = 4 V). Each panel cor- 
respond to the following experimental (/ exp ) and numerical 
(/num) driving frequencies in kHz: (a) (/exp, /num) = (275,268), 
(b) (/cx P ,/nu m )=(276,269), (c) (/ex P , /num) =(280,273), (d) 
(/oxp,/num) = (248,244), (e) (/ oxp , /„ um ) = (254,247) and (f) 

(/oxp,/num) = (249,245). 



than the 3-peaked one. The reason is that in Fig. [T27 f), 
the 3-peak solution was obtained at higher frequencies 
and then adiabatically extended to lower ones, whereas 
in Fig. lE27 e) the 2-peak solution was obtained starting 
from the 1-peak region. 

We show the eventual location of peaks in the breather 
pattern (i.e. after the driver has been on for a long time). 
However, it is important to mention that the exact loca- 
tion where the peaks eventually settle is sensitive to slight 
impurities in the lattice. We have noticed that when 
we turn on the voltage source, at first we can observe a 
more sinusoidal pattern (corresponding to the most mod- 
ulationally unstable fc-value), but as the pattern reaches 
higher energy and becomes more nonlinear, the peaks 
may shift and adjust themselves in the lattice. As it can 
be observed, peaks are not perfectly equispaced in the 
lattice. This is obviously due to the inhomogeneities and 
noise present in the experiment. 

For the numerical results depicted in Fig. [T^] we used a 
set of initials conditions based on the experimental data 



and determined the stationary state by letting the nu- 
merical profiles to settle to a steady state. For the cases 
corresponding to Vd = 2 V and Vd = 4 V, adding a small 
frequency offset A/ rj 4 — 7 kHz, we observe, in general, 
a good agreement between numerics and experimental 
data. The mismatch between experiments and theory, in 
particular the intersite distance peaks, can be attributed 
to the above mentioned factors. Furthermore, to repro- 
duce precisely the experimental peak voltage is extremely 
difficult because it corresponds to the voltage at reso- 
nance and, therefore, even very small parameter changes 
can create large differences in the maximum amplitudes. 
Nevertheless, the quantitative agreement appears fairly 
good, especially for the Vd = 4 V case. 



IV. CONCLUSIONS 

In this paper we have formulated a prototypical model 
that is able to describe the formation of nonlinear in- 
trinsic localized modes (or discrete breathers) in an ex- 
perimental electric line. This has been derived based on 
a combination of the fairly accurate characterization of 
a single element within the lattice (including its nonlin- 
ear resonance curves and hysteretic behavior) and fun- 
damental circuit theory in order to properly couple the 
elements. Comparison between theory and experiments 
shows very good qualitative and even good quantitative 
agreement between the two. We characterized the re- 
gions of existence and stability of n-peaked breathers for 
n = 1,2,3 and illustrated how transitions of the coher- 
ent waveforms of one kind to those of another kind take 
place, rationalizing them on the basis of stability prop- 
erties and their corresponding Floquet spectra. We also 
showed that the precise number of peaks and their loca- 
tion in the lattice is fairly sensitive to initial conditions, a 
feature also generally observed in the experiments where 
the potential for states with different numbers of peaks 
similarly manifests. 

Naturally, many directions of potential future research 
stem from the fundamental modeling and computation 
basis explored in the present manuscript. On the one 
hand, it would be very interesting to attempt to under- 
stand the stability properties of the different breather 
states from a more mathematical perspective, although 
this may admittedly prove a fairly difficult task. On 
the other hand, from the modeling and computation per- 
spective in conjunction with experimental progress, the 
present work paves the way for potentially augmenting 
these systems into higher dimensional setups and at- 
tempting to realize discrete soliton as well as more com- 
plex discrete vortex states therein 0, [l8| . Such studies 
will be deferred to future publications. 
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